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ABSTRACT 

UV irradiation of simple ices is proposed to efficiently produce complex organic species during star- 
and planet-formation. Through a series of laboratory experiments, we investigate the effects of the 
H2O concentration, the dominant ice constituent in space, on the photochemistry of more volatile 
species, especially CH4, in ice mixtures. In the experiments, thin (^^40 ML) ice mixtures, kept at 20- 
60 K, are irradiated under ultra- high vacuum conditions with a broad-band UV hydrogen discharge 
lamp. Photodcstruction cross sections of volatile species (CH4 and NH3) and production efficiencies 
of new species (CaHg, C2H4, CO, H2CO, CH3OH, CH3CHO and CH3CH2OH) in water-containing 
ice mixtures are determined using reflection-absorption infrared spectroscopy during irradiation and 
during a subsequent slow warm-up. The four major effects of increasing the H2O concentration are 
1) an increase of the destruction efficiency of the volatile mixture constituent by up to an order of 
magnitude due to a reduction of back reactions following photodissociation, 2) a shift to products rich 
in oxygen e.g. CH3OH and H2CO, 3) trapping of up to a factor of five more of the formed radicals in the 
ice and 4) a disproportional increase in the diffusion barrier for the OH radical compared to the CH3 
and HCO radicals. The radical diffusion temperature dependencies are consistent with calculated 
H20-radical bond strengths. All the listed effects are potentially important for the production of 
complex organics in H2 0-rich icy grain mantles around protostars and should thus be taken into 
account when modeling ice chemistry. 

Subject headings: astrochemistry; methods: laboratory; ISM: molecules; circumstellar matter; molec- 
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ular data; molecular processes 

1. INTRODUCTION 



Observations of complex molecules toward star-formiuK 



pre-biotic formation mechanisms (e.g. 



Belloche et al 



2009; 



regions and in comets demonstrate the existence of efficient 



Crovisier et al 



was suggested as a path to chemical complexity more than three decades ago 



2004 ). Photochemistry in icy grain mantles 



Greenberg &: Hone 



19741 ). Despite this 



history and a recent surge in interest, most mechanisms that produce complex molecules during UV an d ion irradiation 



of simple ices remain poo rly understood. The photochemistry has been quantified for some pure ices (jGerakines et al 



1996; 



Oberg et al 



2009a|) . In astrophysical settings the chemistry is complicated by the observation that most ices are 
mixed. Predicting the outcome of ice irradiation in space thus requires a quantitative understanding of how different 
mixture constituents affect the ice photochemistry. Because of its prominence and its known effects on the ice binding 
environment, this study focuses on the effect of H2O ice during UV irradiation of binary ice mixtures. These mixtures 
are not proposed to perfectly mimic the multi-component ice mixtures found in space, but are used to probe the 
fundamental principles of ice photochemistry; principles that can then be applied to more complicated ice systems 
through a combination of modeling and further experiments. 
Ices are common during star formation, with II2O ice reaching abundances of 10^* with respect to II2 (e.g. 



Boogert et al 



20081 ). Simple ices - H2O, CO, CO2, NH3, CII4 and CII3OH - form in molecular cloud (cores) through 



direct freeze-out and through hydrogenation and oxygenation of atoms - O, C and N - and of unsat urated molecules 



such as CO on grain surfaces (ITielens fc Hagen 



ice spectral features and ice maps 



Bergin et al 



1982 



Watanabe fc Kouchi 



2002; 



Pontoppidan et al 



2002' 



2003 : 



loppolo et al 



Pontoppidan 



2008f) . From observed 



20061 ). the formation of 



ices is sequential, starting with H2O and CO2 (and probably CH4 and NII3 as well). Deeper into the cloud core or 
later during the core contraction CO freezes out catastrophically, resulting in a s econd ice phase dom inated by CO 



Oberg et al 



(l2009al) . and of CH4 



and later CII3OH. The photochemistry of CHsOHiCO ice mixtures, investigated in 
NII3 and CO2 in H2O ice mixtures are then prime targets for laboratory experiments. 
Irradiation of CH3OH and CH30H:C0 ices result in the formation most complex C,H ,0-bearing species, such as 



2009al ). The same molecules may 



HCOOCH3, CH3CHO and C2H5OH, found in protosteUar hot cores (jOberg et al. 
also form from irradiation of CH4:H20 ice mixtures, however, and the relative contribution from the two ice phases 
to the observed complex molecule abundances are unknown. Understanding the chemistry in the H20-rich ice is 
furthermore required to predict the formation rates of the pre-biotically important N-bearing complex molecules, since 
the main source of N in the ice, NH3/NH^, probably forms mixed with H2O (Bottinelli et al. ApJ in press). 
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The response of ice mixtures to irradiation also governs the amount of e.g. CH4 that remains for thermal desorption 
close to the protostar. T hermal desorption of CH4 ice is thought to result in a complex warm-carbon-chain-chemistry 



around some protostars (jSakai et al. 



20081) . One of the aims here is to investigate how the CH4 ice photodestruction 



rate depends on H2 0-ice concentration and ice temperature. 
Irradiation of H20-rich ices was first investigated in the 60s and 70s, (e.g. 



Hagen et al 



19791) . Since then there 



has been a handful of studies that have included the effects of H2O on the overall photochemistry and two dedicated 
studies on the effect of different H2O concentrations on the CH4 chem istry during prot o n bombardment a t 10 K 



(Moore fc: Hudson 



19981) and on the CH3OH UV photochemistry at 3 K (jKrim et al 



2009) 



Moore fc Hudson! (|l998l ) 



found an increasing production of CH3O H, CH3CHO and CH3CH2OH and a decreasing production of C2Hg with 



H2O concentration in H20;CH4 mixtures. 



Krim et al 



(j2009f) found an almost constant CH3OH conversion into CO, 



CO2 and H2CO for a pure and a H20:CH30H 1:1 ice, but a factor of 3-7 conversion increase in a H20:CH30H 10:1 
mixture. In both studies all changes were explained by the production of increasing amounts of OH radicals with H2O 
concentration and subsequent radical-radical reactions or H-abstraction by OH. 
H2O may however affect the photochemistry in more ways than increasing the number count of OH radical s in the 



ice. The binding energies in H2 0-rich ices are typically different compared to pure ices ([Collings et al. 



20031 ) ■ which 



may affect radical diffusion. R adicals may also bec ome physically trapped in H2O ice at low temperatures as is often 



observed for volatile molecules (jCollings et al 



120041) . These radical-H20 interactions have probably different strengths 
for different radicals, which may drive the chemistry in otherwise unexpected directions. 

We investigate the relative importances of these potential effects of H2O on the photochemistry of CH4, NH3 and 
CO2 ices. The focus is on quantitative comparisons between pure and binary H20:CH4 ices at different concentrations 
and temperatures, both during irradiation and during the subsequent warm-up. 

2. EXPERIMENTS 



20051 ) under ultra-high vacuum conditions 



The experiments are carried out on the set-up CRYOPAD (jOberg et al. 
(^10~^-10~^*' mbar). Ices are deposited diffusively at 18 K by introducing a gas (mixture) in the vacuum chamber 
along the surface normal of a gold substrate, which is temperature controlled down to 18 K with a 2 K uncertainty. 
10-20 mbar gas mixtures are prepared in a separate glass manifold with a base pressure of lO^'' mbar. The CH4, NH3 
and "'^'^002 (^"^002 was used to minimize overlap with gas phase CO2 features from outside of the vacuum chamber) 
gases have a minimum purity of 99.9% (Indugas). Samples containing H2O are prepared from the vapor pressure of 
de-ionized H2O, further purified through freeze-thawing. 

The original ice mixture as well as changes in the ice composition induced by UV irradiation are quantified through in- 
frared spectroscopy in refiection-absorption mode (RAIRS). The relative RAIRS band strengths are consis tent with rel 



ative transmission band strengths in the investigated ice thickness regime and thus certain within 20-30% ( Oberg et al 



Oberg et al. 



2009g). Absolute band strengths have a 50% uncertainty, but this does not affect the quantification of the chemistry, 
where yields are calculated in fractions of the original ice. 
The ices are irradiated by a hydrogen-discharge UV lamp, peaking at Ly-a but extending between 6 and 11.5 eV 



(Mufioz Caro &: Schutte 



20031 ). Ah ices are irradiated with a UV flux of -1.1(±G.5) x lO^^ cm^^ fo^. g hours, 
resulting in a total flucnce of ^2.3x10^^ cm^'^. This is comparable to the fluence an ice in a cloud core is exposed to 



during 10^ years because o f cosmic-ray induce d UV photons at a flux of lO"' cm ^ s ^ ( Shen et al 



calibration is described in ( Oberg et al 



2009c ) 



20041) ■ The lamp 



Table [T] lists the photochemistry experiments in terms of their mixture composition, ice temperature during UV 
irradiation and total thickness. Following irradiation at 20-100 K, all experiments are heated by 1 K min^^ to 150- 
250 K, while acquiring RAIRS every 10 min. The ice thicknesses range between ~15-54 monolayers (ML), but most 
experiments are carried out with '-~^40 ML ices. This ice thicknesses are similar to what is expected in the dense and 
cold stages of star formation and they are also in the linear regime for RAIRS. The ices are also thick enough that 
photod esorption wiU not significantly affect the bulk of the ice, based on previously measured UV photodesorption 



yields (jOberg et al 



2009cllbl ) 



3. RESULTS 



The following five sub-sections focus on 1) the photodestruction cross sections of CH4, NH3 and CO2, 2) the identi- 
fication of the photoproducts, 3) the H20:CH4 photochemistry as a function of H2O concentration during irradiation 
and (4) during warm-up, and 5) the effects of temperature on the H20:CH4 UV photochemistry in dilute mixtures. 



3.1. Photodestruction cross sections in pure and mixed ices 

The photodestruction rate of a species determines the production rate of radicals in the ice. Its maximum value 
is the photodissociation rate, as measured in the gas-phase. Previous ice experiments have however sh own that the 



meas ured ice photodestruction rates are substantially lower than the gas-phase photodissociation rates (jCottin et al 



20031) . probably because of fast back-reactions between the dissociation fragments. Figure [T] shows the ^'^C02, CII4 
and NH3 spectral features before and after a UV fluence of 2.3 x 10^^ cm~^ in pure ices and Il20:X ^5:1 ice mixtures 
at 20-30 K. During irradiation, more CII4 and NH3 are destroyed in the mixtures than in the pure ices. 

This effect is shown quantitatively in Fig. [51 where the loss of CH4 and NH3 in the different ices is plotted as a 
function of UV fluence (overlap with gas-phase CO2 lines prevents a similar analysis for CO2 ice). The photodestruction 
cross sections are calculated from the first 4 x 10^^ photons cm~'^, where the curves are still approximately linear, 
and the results are listed in Table [1] - the listed uncertainties are absolute, the fit uncertainties are 10-20%. The 
photodestruction cross sections for pure CH4 and NH3 are 0.5 and 1.4x10"^^ c m^, respectively. Bo th are an order 



of magnitude lower than the measured gas phase photodissociation cross sections (jvan Dishoeck 



19881) . They are still 



I 

higher than previously reported values for 1000 ML ices (jCottin et al 



20031) . suggesting that those experiments may 



have suffered from optical depth effects. 

The photodestruction cross sections are up to an order of magnitude higher in the H2O ice mixtures, close to the gas 
phase photodissociation values. In the case of CH4, even the small amount of H2O in the H20:CH4 1:3 ice mixture 
results in a factor of five higher destruction cross section, while a H20-dominated ice (5:1) is required to increase the 
NH3 cross section significantly. 

The increasing UV photodestruction with H2O concentration suggest that H2O molecules surrounding the photodis- 
sociated volatiles may trap some of the radicals and thus inhibit reformation of CH4 and NH3. The reason behind 
the differences between different species is difficult to assess without further modeling, but the effect seems larger, 
the lighter the volatile dissociation radicals (CH3+H vs. NH2+H vs. CO+0), assuming that the lack of an effect of 
H2O on the final CO2 dissociation fraction corresponds to a constant CO2 photodissociation cross section with H2O 
concentration. 

Temperature does not affect the CH4 photodestruction cross section between 20 and 60 K for dilute (^^5:1) H20:CH4 
ice mixtures (Fig. [2};) . The effect on more CH4-rich ice mixtures cannot be tested because of thermal outgassing of 
CH4 above 30 K when in high concentrations. In comparison, recent experiments on pure CH3OH ice photodissociation 
resulted in an increasing photodes truction (then termed 'effective photodissociation') cross section with temperature 



by up to 50% (jOberg et al 



l2009al) . In the CH3OH study this was explained by an increased escape probability of the 
dissociation fragments at higher temperatures, thus lowering the back-reaction rate. A similar effect was expected for 
CII4. Its absence suggests that thermal diffusion is fast already at 20 K for the volatile CH4 dissociation fragments 
compared to the main CII3OH photodissociation product, CH2OII. 

3.2. Product identification 

Quantification of the ice photochemistry requires secure identifications of the main reaction products. This was not 
possible for the UV-irradiated H20:NIl3 and H20:C02 ices, because the spectral features of expected main products - 
NH2OH, HCOOH, II2CO3 - are blended with bands of II2O, NH3 or other photoproducts. The expected photoproducts 
in Il20:CH4 ices are shown in Fig. [3] Figure 2] shows that most of the spectral features appearing upon irradiation 
of Il20:CH4 ice mixtures can be assigned to th e predicted main produc ts: C2H6, II2CO, CO (not shown), CII3OH, 



CII3CH2OII and CH3CHO, in agreement with iMoore fc Hudson 



(|l998f ). A few bands remain unassigned and are 



probably due to more complex molecules and to radicals. The figure also shows that two of the products, CII3OII 
and CII3CII2OII, have bands that shift significantly between the pure ice and the Il20-mixture. This introduces some 
uncertainty in the assignment of CH3CII2OH - it is only possible to securely separate the relative contributions of 
CH3OH and CH3CH2OII in a subset of the experiments. 
In general the spectral bands used for identification and quantitative analysis have minimum overlap with spectral 
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2009a|) . These band positions 



features of other detected species as weh as of more complex molecules (jOberg et al. 
and strengths are listed in Table [2l When deriving ice abundances, the region of the investigated bands is typically 
fitted with several Gaussians to account for band overlaps. This approach was chosen instead of fitting ice spectra 
because of the varying spectral shapes in different ice environments. The same band strengths were adopted for all ice 
mixtures, since the exact band strengths are unknown for most of the used ice mixtures. 

3.3. H2 O concentration effects at 20 K 

Figure [5] shows that irradiation of pure CH4 ice results in the production of C2H6, C2H4 and larger hydrocarbons, as 
expected - no C2H2 was observed and in general the products are hydrogen-rich. Adding H2O to the CH4 ice results 
in the appearance of CO, H2CO, CH3OH, CH3CHO and CH3CH2OH spectral features. The same bands are visible 
in all H2O ice mixtures, but the relative contributions of different bands change with concentration. 

The effect of H2O concentration on the final product abundances, with respect to CH4, is quantified in Fig. [51 
The identified products can be sorted into three groups a) hydrocarbons which have increasing abundances with CH4 
concentration, b) small organics like CH3OH that form from a CH4 and a H2O fragment, and c) larger organics like 
CH3CH2OH that form from two CH4 fragments and one H2O fragment. The molecules in each group are then expected 
to form as: 



A^CH 

A^CHsOH 



^''^"^ OC [CH4], (1) 



AcH4 
-^CHaCHaOH 



OC 



(100 - [CH4]) and (2) 



OC (100 - [CH4])[CH4]. (3) 



A^CH4 

Qualitatively the abundance patterns of these three molecules agree with the curve predictions. The difference in 
e.g. the peak position of the curve in Fig. |BJ: indicates that the effects of diffusion and different photodestruction 
cross sections for different mixtures cannot be neglected when making quantified predictions, however. In addition, the 
production dependencies of H2CO, CO and CH3CHO are expected to be complicated by multiple formation pathways, 
all involving multiple dissociation events as shown in Fig. |31 

Four of the detected photoproducts, C2H6, CH3OH, H2CO and CO, form abundantly enough to quantify their 
production as a function of fluence during UV irradiation (Fig. [7]). C2Hg and CH3OH form immediately upon 
irradiation in all experiments. The H2CO production is delayed in the H20-dominated mixture and CO production 
does not begin before 5 x 10^^ photons cm^^ in any mixture as expected for a reaction pathway through multiple 
photodissociation events. 

The initial photoproduction rates of C2H6, CH3OH and H2CO are calculated from the abundance growth during 
the first 4 x 10^'^ photons cm when the growth is still roughly linear. Table [3] lists the resulting formation 'cross 



7 

sections'. These numbers do not represent physical photoproduction cross sections, since diffusion is required for their 
formation, it is however a convenient way of parameterizing their formation. The total formation cross section of 
complex molecules, in number of product molecules with respect to CH4, increases somewhat with H2O concentration. 
When taking into account that some products require multiple CHx fragments, an approximately constant fraction of 
the original CH4 is incorporated into more complex molecules at all concentrations. Furthermore, the CH4 photode- 
struction cross section increases with an order of magnitude, between pure CH4 and the 4:1 mixture and the formation 
efficiencies normalized to the photodestruction cross sections thus decrease significantly with H2O concentration. The 
presence of H2O then slows down the chemistry of the radicals in the ice, even though this effect is 'compensated for' 
under laboratory time scales by the even more efficient slow-down of back reactions into CH4. 

3.4. Warm up of irradiated ices 

During warm-up of the UV- irradiated ices, formation of new species together with sequential desorption change the ice 
mixture with temperature (Fig. [5]). In the H20-poor 1:3 mixture CH3OH is the only molecule that forms perceptively 
during warm-up, in the 2:1 mixture CH3CHO is alone to form until reformation of the H20-ice network around 
130 K. The different formation patterns in the H2 0-poor and H2 0-rich ices suggest different diffusion environments. 
Specifically, OH diffusion seems hindered in the H2 0-rich ice, while HCO can still diffuse and react with CH3 to form 
CH3CHO. All species are partially trapped in the H2O ice, so that even C2H6 desorption is only complete around the 
H2O desorption temperature of ^150 K 

3.5. Ice temperature effects during irradiation at 20-60 K 

Increasing the ice temperature should speed up the thermal diffusion of radicals in the ice and thus the ice chemistry. 
Figure |9] shows that the most notable changes in the irradiated ice spectra at 20, 40 and 60 K are instead due to 
outgassing or a decreased production of volatile species, such as H2CO and CO, at the higher temperatures. The 
decreased production of these secondary and tertiary dissociation products can be understood if the intermediate 
radicals, e.g. HCO diffuse fast enough at higher temperatures to react with another radical before absorbing a second 
photon and dissociating to CO. This scenario is supported by an earlier onset in the CH3CHO production at 60 K 
compared to 20 K (Fig. [TU)) . and also consistent with the warm-up results that reveal HCO diffusion above 30 K. 

In contrast, there is no evidence for an increased formation efficiency of CH3OH and CH3CH2OH or for the formation 
of more complex molecules such as (CH20H)2. Rather the CH3OH production cross section decreases from 5.6 to 
3.3 10^^® cm^ / NcH4(0) between 20 and 60 K. The lack of a positive temperature effect on the CH3OH production, 
illustrated in Fig. [TUl suggests that CH3OH and CH3CH2OH formation depends mainly on CH3 diffusion that is 
efficient already at 20 K; OH and CH2OH diffusion must be too slow, even at 60 K, in a H20-dominated ice to affect 
their formation. In addition the decreasing cross section with temperature is probably due to thermal desorption of 
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some of the produced CH3 at 60 K before it has time to react. 



4. CALCULATIONS OF BINDING ENERGIES 



To aid the interpretation of the experiments, we calculated the binding energies of radical-H20 complexes using 
the CCS D(T) (coupled cluster with single, double, and perturbative triple excitations) method with the aug-cc-pVTZ 



basis set ( Kendall et al 



counterpoise method (jBovs fc Bernardi 



1992). In order to c orrect for the Basis Set Superposition Error (BSSE) the Boys-Bernardi 



1970r) was used. The internal geometries of the interacting species were kept 
constant while a large number of different relative orientations and intermolecular distances were scanned. The thus 
located lowest energy minima were used to calculate the binding energies. The reported energies are electronic energies 
relative to the infinitely separated interaction partners. No account of zero -point energy was made. All calculations 



were performed using the Gaussian 03 program package (jFrisch et al. 



2004) 



The calculations were carried out for the predicted three most common radicals in the ice, OH, HCO and CH3. The 
resuhing binding energies are 0.23 eV (2700 K) for H2O-OH, 0.11 eV (1300 K) for HaOtHCO and 0.07 eV (800 K) for 



H2O-CH3. The H2O-OH result is in good agreement with other recent high-level calculations (jSoloveichik et al 



20101) 



5. DISCUSSION 
5.1. The influence of H2O 

Under laboratory conditions the most important effect of increasing the H2O concentration is to drive the chemistry 
away fro m Cxtly chains and towa rd products with a higher and higher concentration of OH. The same effect was 



noted by iMoore fc HudsonI ()1998[ ) during proton bombardment of H20:CH4 ice mixtures. This can be understood 
from simple rate equations where the production rates of various molecules are proportional to the production rates 
of CH3/CH2 and OH in the ice, i.e. the photodestruction cross sections of CH4 and H2O. The absolute values of these 
reaction rates also depend on the diffusion rates of radicals. If the destruction and diffusion rates had been constant, 
the relations shown in Fig. [5] could have been used to derive the relative diffusion barriers directly. 

H2O does, however, change the CH4 photodestruction cross section. Since this increased destruction efficiency is not 
accompanied by an equally large increase in production efficiencies of new molecules, H2O must slow down the overall 
diffusion of radicals in the ice to the point of trapping them. H is not expected to be trapped at 20 K, however, and 
the existence of larger radicals in the ice indicates that many of the H atoms formed during dissociation either escapes 
from the ice or react to form H2. A significant amount of the produced heavier radicals, especially CH3, must also 
remain mobile in or on top of the ice to account for the production of any complex molecules during irradiation at 20 K 
in the H2 0-rich ices. Without detailed modeling it is not possible to estimate the relative importances of surface and 
bulk diffusion for the chemistry. It is however unlikely that all radicals are trapped in the bulk of the ice at 20 K, since 
even in the most dilute mixtures ~20% of the dissociated CH4 molecules recombine to form more complex molecules. 



This corresponds to at least 8 ML of chemically active ice in a 40 ML thick ice and the ice would thus have to be very 
porous to explain all chemistry with diffusion on external and internal (in large pores) ice surfaces. It would also be 
strange if diffusion in H2O ices is fundamentally different comp ared to diffusion in CH3OH ices, where bulk diffusion 



is definitely required to explain the ice ph otochemistry result s ( Oberg et al 



of the two ices during segregation studies ([Ehrenfreund et al. 



199S; 



2Q09a ) , cons idering the similar behavior 



Bernstein et al 



20051) 



Trapping of some of the radicals in the ice bulk is however needed to explain that up to 80% of the dissociated 
CH4 is not converted into more complex molecules in the H20:CH4 5:1 ices. This is con sistent with the commonly 



20041 ) . It may also explai n some 



Krim et al 



(|2009l ) would 



observed trapping of molecules such as CO and CO2 in Il2 0-rich ices (jCollings et al. 
earlier results; trapping of CII3OII photofragments in the H20:CIl30II 10:1 ice mixture in 
prevent the radicals reacting with each other and thus favor dissociation into smaller and smaller fragments (CH2OH, 
II2CO, HCO, CO), explaining their result of high H2CO and CO yields in H20-rich ices without invoking efficient H 
abstraction by OH radicals. 

When the ice is heated some of the bound radicals overcome their entrapment and react to form new species. In 
all H20:CIl4 ice mixtures the most abundant radicals are expected to be CII3 (see below) and OH, the two major 
photoproducts of CH4 and H2O. CH3OH is thus always expected to form more abundantly during warm-up of the 
irradiated ices than e.g. CH3CHO, consistent with the results in Fig. [7] for the H20:CH4 1:3 experiment. In the 
2:1 mixture almost no CH3OH forms initially during warm-up, however. Instead a significant amount of CH3CHO 
forms, followed later by CH3OH formation during the H2O ice reformation temperature of ^^120-130 K. An increasing 
H2O concentration thus increases the OH diffusion barrier relative to the HCO diffusion barrier. This is physically 
reasonable since OH and H2O are expected to bond stronger than HCO and H2O. The importance of the relative 
diffusion barriers of OH and HCO to drive the chemistry demonstrates that radical-radical reactions in ice is only 
efficient when both radicals (i.e. CH3+OH or CH3+HCO) taking part in the reaction arc mobile. 

5.2. Toward quantifying the H20:CHi photochemistry 

CH4 is known to dissociate into a number of different fragments during UV photolysis in the gas phase. The 
dissociation cr oss section is wavelen gth dependent, e.g. almost no CH2 is observed to form from dissociation with 



Ly-g photons (IMordaunt et al 



(jSlanger fc Black 



19931 ) even though it is expected to be an important channel at lower photon energies 



19821 ). By comparing the initial C2Hg and C2H4 formation efficiencies in pure UV-irradiated CH4 



ice, it is possible to constrain the CH4 photodissociation branching ratio in ices irradiated with a hydrogen lamp output 
consisting of both Ly-a photons and broad-band emission. 

In the pure CH4 ice the initial growth of both C2He and C2H4 is assumed to be caused by CH3-I-CH3 and CH2+CH2 
radical reactions, ignoring second generation photodissociation products. Further assuming that the diffusion barriers 
of CH3 and CH2 are comparable, the CH4 branching ratio can be constrained from the initial C2H6 and C2H4 



10 



production rate, which is ~9:1; 



Gerakines et al 



Oberg et al. 



(|l996f) find a comparable C2H6/C2H4 product ratio (8:1). Since each 
reaction requires two radicals, the dissociation branching ratio is inferred to be CH3:CH2 ^3:1. This is consistent with 
a chemistry driven largely by Ly-a photons. 

The diffusion barriers are more difficult to constrain directly from the experiment. The calculated bond strengths 
of 800, 1300 and 2700 K for the H2O-CH3, -HCO and -OH complexes are not identical to ice diffusion barriers, but 
the relative values should agree with the experimental results. Indeed, the data show that OH has a higher diffusion 
barrier than HCO in H20-dominated ices (efficient thermal diffusion at 30 and 120 K respectively) and that CH3 must 
be diffusing already at 20 K to explain any complex molecule production at this temperature. The calculated bond 
strengths are thus a good starting point for models, where the quantified formation rates of molecules can be used to 
extract diffusion barriers, similarly to what is currently being done for the CH3OH photochemistry (Garrod & Oberg 
in prep.). Already both calculations and experiments suggest that the temperature window between diffusion and 
desorption in H20-dominated ices is quite small for the investigated radicals, which may explain the low production 
yield of complex organics during warm-up (a few % of all dissociated CH4) observed here compared to the efficient 



production of complex molecules observed during warm-up of irradiated CH30H-ices (|Oberg et al 



2009a ) 



The calculations also show that predicting radical d iffusion barriers from the bond strengths of the parent molecule 



to H2O can be quite reasonable (jGarrod et al 



0.2, 0.07-0.13 and 0.03 eV, respectively (e.g. 



20081): the H9O-H9O H2O-H2CO and H2O-CH4 bond strengths are 



Szczesniak et al 



1993 : 



Bene 



1973I ). which is very similar to the bond 



strengths calculated for the H2 0-radical complex interactions. There thus seems to be a direct dependence between 
the H2 0-molecule and H2 0-radical interactions. 

5.3. Astrophysical implications 

As shown previously for pure and mixed CH3OH ices, UV irradiation of H20:CH4 ices readily results in the pro- 
duction of complex molecules and this production can be quantified. After a UV fluence corresponding to ^ 10^ years 
in a cloud core, up to 50% of the original CH4 and 25% of the NH3 ice have dissociated into radicals that can react 
into more complex species. At low temperatures most of these radicals are trapped in the ice, i.e. only ^10% of the 
original CH4 ice can be converted into complex molecules during UV irradiation at low temperatures. As the protostar 
turns on and heats the ice many of the trapped radicals will become mobile resulting in a second equally important 
formation step of complex molecules - this step will probably be more important in space than in the lab because of 
the lower heating rate in astrophysical environments. Taking a typical CH4 ice abundance of ^4% with respect to H2O 



ice and a H2O ice abundance of ^^10 with respect to H2 toward protostars (jOberg et al 



20081) . this corresponds to 



a potential production of complex molecules of 0.04 x 10 ^ x 0.1 = 4 x 10 ^ with respect to H 2. This is high enoug h 



to contribute significantly to typical hot core abundances of complex molecules ~ 10 ^ — 10 ^ (|Bisschop et al 



2007 ) 



Several of the CH3OH and H20:CH4 photochemistry products are the same, e.g. CH3CH2OH and CH3CHO, and 
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the relative importance of the two production pathways during star formation wih depend on the relative abundances 
of CH3OH and CH4 ice on the grains as well as what other ice constituents are mixed with H2O. Both the H2O- and 
the CO/CHsOH-rich ice phases must clearly be modeled to accurately predict the total formation rates of C- and 
0-containing complex molecules in space. 

In addition it is crucial to understand the effects of H2O quantitatively to constrain when and where N-containing 
complex molecules form - their main formation pathway is probably from NH3 in a H2 0-rich ice. While the production 
of radicals seems to be the most important effect of H2O on the ice photochemistry under laboratory conditions, this is 
not necessarily the case for the time scales and ice compositions present at star and planet formation. Rather detailed 
modeling is required that extracts the microscopic properties of the H20-rich ice chemistry from the laboratory 
experiments, such as the composition-dependent diffusion barriers, and applies these results to astrochemical models. 
Only when a simple system, such as H20:CH4, is understood at this detailed level can we expect to accurately model 
the more elusive nitrogen chemistry in H2 0-rich environments and thus provide predictions of e.g. the prebiotic amino 
acid production. 

Estimating ice life times is easier than predicting the chemical evolution, though the composition dependence of 
CH4 photodestruction introduces some uncertainty. Table 2] list the photodestruction time scales for 20-30 K CH4 
ice in pure, 2:1 and 4:1 H20:CH4 ices subject to a weak UV field induced by cosmic rays, th e interstellar radiation 



19951 ). Outside of protected 



field and the 1000 times higher UV fiux inferred toward the LI 527 outflow (jSpaans et al. 
cloud cores the CH4 ice life time is thus short compared to most other time scales, especially in H2 0-rich ice mixtures 
relevant for protostars. For large quantities o f CH4 ice to surv ive to thermally desorb around protostars and drive 



warm-carbon-chain-chemistry as suggested by (|Sakai et al. 



20081) then requires that the ices are well protected up until 



that point. This may be a reason for the low number of observations of such a chemistry during the protostellar stage. 

6. CONCLUSIONS 

1. The photodestruction cross section of NH3 and CH4 in H2O ice mixtures increases with H2O concentration by 
up to a factor of 10 with implications for the CH4 ice lifetime in protostellar envelopes. Simultaneously the 
production efficiency of stable molecules in the H20:CH4 ice increases by a factor of two. This is explained by 
trapping of radicals in the H20-ice matrix, which prevents back reactions between e.g. CH3 and H. 

2. The H20:CH4 photochemistry produces C2H6, CH3OH and CH3CH2OH as a function of CH4 concentration, in 
a way that can be directly related to the number of OH and CH2/CH3 groups in each product. The production 
of CH3CHO, H2CO and CO has a fluence delay consistent with formation through more steps compared to e.g. 
CH3OH. 



3. On laboratory time scales increasing the ice temperature from 20 to 60 K has a limited effect on the chemistry. 
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probably because HCO is the only important reactant that becomes mobile in this temperature interval in 
H20-ciominated ice mixtures. 

4. The CH4 photodissociation branching ratio into CH3 and CH2 is >3:1, as derived from the relative production 
efficiencies of C2H6 and C2H4. 

5. OH diffusion is fast in the H20:CH4 1:3 mixture, but not present in the H20-dominated ice mixtures below 
100 K. HCO diffusion is possible above 30 K in the H20-rich ice mixtures and CH3 diffusion seems efficient 
at 20 K in all mixtures, consistent with calculated H20-radical complex bond strengths. OH diffusion is thus 
more severely slowed down with H2O concentration than the other radicals, which counteracts the higher OH 

production rate in H20-richer ices. 

6. Because H2O is both a source of OH radicals and preferentially trap OH compared to more volatile radicals, 
predicting the photoproduction branching ratio of different complex molecules in space requires microscopic 
modeling of the H20:CH4 experiments rather than direct comparison with experiments. 

This work has benefitted from discussions with Herma Cuppen and Robin Garrod. Support for KIO is provided 
by NASA through Hubble Fellowship grant awarded by the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, Inc., for NASA. Astrochemistry in Leiden is supported 
by a SPINOZA grant of the Netherlands Organization for Scientific Research (NWO) and NOVA. The work by SA 
was partly funded by a NWO-CW Top grant. Some of the electronic structure calculations were performed at C3SE 
(Chalmers Centre for Computational Science and Engineering) computing resources. 
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TABLE 1 
The photochemistry experiments. 



Ice 


H20:x 


Temp. (K) 


Thick. (ML) 


cross section (10 


CH4 


0:1 


30 


47 


0.5[0.3] 


H20:CH4 


1:3 


20 


42 


2.8[1.4] 


H20:CH4 


2:1 


20 


37 


2.9[1.5] 


H20:CH4 


4:1 


20 


37(29:8) 


4.9[2.5] 


H20:CH4 


4:1 


40 


38(30:8) 


4.6[2.3] 


H20:CH4 


5:1 


60 


38(32:6) 


3.9[2.0] 


NH3 


0:1 


30 


51 


1.4[0.7] 


H20:NH3 


1:1 


20 


54 


1.5[0.8] 


H20:NH3 


4:1 


20 


43 


5.0[2.5] 


CO2 


0:1 


30 


15 




H20:C02 


6:1 


20 


35 
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TABLE 2 

Ice infrared spectral feati'res i'sed for quantification. 



Species 


Band (cm ^) 


Band strength'' (cm ^) 


Reference 




CH4 


1300 


6.1x10"^* 


Moore & Hudson (1998) 




NH3 


1070 


1.7x10"" 


D'Hendecourt & AUamandola ( 


1986) 


CO2 


2343 


7.6x10"" 


Gerakines et al. (1995) 




H2O 


1670 


1.2x10"" 


Gerakines et al. (1995) 




C2H6 


2976 


1.1x10"" 


Moore & Hudson (1998)'' 






821 


1.9x10"^* 


Pearl et al. (1991) 




C2H4 


1436 


2.9x10"^* 


Moore & Hudson (1998)'' 




H2CO 


1500 


3.9x10"^* 


Schutte et al. ( 1993) 




CH3OH 


1026 


2.8x10"" 


D'Hendecourt & AUamandola ( 


1986) 


CH3CH2OH 


1044 


7.3x10"^* 


Moore & Hudson (1998)'' 




CH3CHO 


1350 


6.1x10"^* 


Moore & Hudson (1998)'' 




CO 


2139 


1.1x10"" 


Gerakines et al. (1995) 





The band strengths are known within ~20-30% when comparing results from different references, ice mixtures and ice 

temperatures. 
'' In a H2O ice matrix. 
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Fig. 5. Difference spectra of CH4 ice and H20:CIf4 ice mixtures after a UV fiuence of 2.3 X 10^'' photons cm ^. The marked peaks 
are a) C2H6, b) CO, c) XCHO, d) H2O, e) H2CO, f) CH3OH, g) C2H6, h) C2H4, i) CH3CHO, j) CH4, k CH3OH, 1) CH3CH2OH, m) 
CH3OH and n) C2H4. 



TABLE .3 

Normalized product 'cross sections' (10""^^ om^ / Nch4(0)) at 20 K for H20:CH4 ices. Fit uncertainties are in brackets. 



Product 


Pure CH4 


1:3 


2:1 


4:1 


C2H6 


0.30[0.10] 


0.52(0.18] 


0.29(0.10] 


<0.20 


C2H4 


0.04[0.01] 








CH3OH 




0.11(0.04] 


0.23(0.08] 


0.56(0.19] 


H2CO 




<0.02 


0.15(0.05] 


0.10(0.04] 


Total" 


0.34 


0.63 


0.67 


<0.86 


Total / Vo"^ 


~100 


41 


33 


^22 



" The total formation cross section of new molecules. 
% of the photodestruction cross section that is accounted for by the total formation cross section taking into account that 

some molecules require multiple CHj, fragments to form. 
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TABLE 4 

CHl If'E LIFE TIMES IX YEARS, 

UV flux (cm-2 s"^) pure CH4 2:1 4:1 

10* (cloud core) 6 x 10® 1 x 10® 6 x 10^ 

10* (cloud edge) 6 x 10^ 1 x 10^ 6 x 10^ 

10" (outflow cavity) 0.6 0.1 0.06 
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Fig. 10.— The quantified production of C2H6, CH3OH, CH3CHO, H2CO and CH3CH2OH in H20:CH4 4-6:1 mixtures at 20 K (stars) 
and 60 K (red diamonds) as a function of UV fluence. The Unes shows the 2-3 data point averages. 



